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ESTIMATION OF FLOOD PROBABILITIES 


Leo R. Beard, A.M., ASCE 


SYNOPSIS 


Data of annual maximum mean-daily flows at 159 long-record river- 
measurement stations are analyzed to determine the characteristics of the 
frequency distributions and the reliability of probability determinations mace 
therefrom. It is shown that the logarithms of the flows are normally dis- 
tributed, which fact makes applicable to flood-flow data a wealth of statistical 
procedures and criteria. An important result of the verification of this fact 
is the ease with which region-wide correlation analyses can be made, through 
which precise hydrologic determinations will be possible. 


INTRODUCTION 


The determination of exceedence probabilities of flood magnitudes is of 
fundamental importance in many design problems. Where adequate records 
are available, the simplest method of estimating such probabilities is by a 
frequency analysis of recorded flood flows. However, the method is presently 
limited in applicability because of the unreliability of extreme-frequency 
determinations. It is the purpose of the following analysis to measure the 
reliability of calculated flood-flow probabilities and to investigate means of 
increasing that reliability. 

In the course of this investigation, it has become apparent that a compre- 
hensive, nation-wide analysis of flood flows will be required in order to es- 
tablish adequate criteria and procedures for flood-frequency determinations. 
The all-important prerequisite to such an analysis is the summarization of 
flood-frequency attributes common to all locations in the form of a charac- 
teristic function of flood-flow distributions. Although many theories have 
been advanced for determining this function, the writer knows of no theo- 
retically complete and adequately documented exposition of a flood-flow dis- 
tribution theory. Information furnished in the following analysis is presented 
with the hope that it will complete the accomplishment of this preliminary 
task insofar as peak flows are concerned. 


Theoretical Distribution of Flood Flows 


By far the most common function that is characteristic of continuous dis- 
tributions of observed values, each value of which is the effect of the identi- 
cal multitude of varying causes, is the normal distribution.! There have 
been many demonstrations that some distributions are fitted by the theoretical 
curve to virtual perfection. For this reason, the normal curve should ordi- 
narily provide the best starting point in a search for the most representative 
function of a continuous variable. 
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In order for a variable to comply with the requirements of a normal dis- 
tribution, it must be continuous, values near the mean must occur more fre- 
quently than extreme values, and it must range from minus infinity to plus 
infinity. The logarithms of flood flows are the values that most logically 
might comply with these requirements. Furthermore, for many years, there 
has been recognized a tendency for the logarithms of annual maximum flood 
flows to be normally distributed.? It is proposed to test this theory by use 
of the wealth of stream flow data collected and published by the U.S. Geologi- 
cal Survey. Inasmuch as characteristics of small samples vary fortuitously 
from those of the parent population, the principal problem to be solved is the 
determination of whether observed variations from the normal are sampling 
variations or are characteristic of the parent distribution (i.e. whether or 
not they can be expected to recur in a future series of events). 


Analysis of Goodness of Fit 
The first step in the analysis of goodness of fit was the selection of 159 


gaging stations, whose records are published regularly in U.S.G.S. Water 
Supply Papers. None of these stations were downstream of substantial con- 


trol or diversion works during the period of record used. With few excep- 


tions, all stations have 25 or more years of record and are located such that 
the center of tributary area is at least 50 miles from the center of tributary 
area for any other station used. Where drainage areas differed by a factor 
of 10 or more, the requirement of 50-mile separation was waived. Stations 
used are listed in Table 1, and locations of tributary-area centers are shown 
in Figure 1. 

Secondly, a tabulation was made of all annual maximum mean-daily flows 
that occurred during the periods of record. Special caution was observed in 
order that values were not excluded because they were too large to record. 

If a gage was destroyed by a flood, an estimate from high-water marks was 
usually available. A few stations where reliable estimates were not available 
were eliminated from the original list, reducing it to the final list of 159 
Stations. 

The third step consisted of tabulating the logarithms of all observed annual 
maximum mean-daily flows, and calculating the mean and standard deviation 
for each station. Standard deviations were computed in accordance with the 
equation: 


3 -/2 (1) 


n-1l 


in which x = deviation of a given logarithm from the mean of the logarithms 
ior the station 
n = number of years of record. 


An example of this computation is given in Table 2. Capital X represents 
the logarithm of a single flow, and lower case x represents the deviation of 
a single logarithm from the mean (M) of the logarithms. The sum of the 
logarithms and the sum of the squares of the logarithms are obtained ina 
single cumulative operation of an automatic calculator. The square of the 
sum divided by the number of events is subtracted from the sum of the 
squares in order to obtain the sum of the squared deviations from the mean. 
This is simply a short-cut operation, which eliminates a large amount of 
tedious calculation. Note that in this case it was necessary to carry 7 places 
in the intermediate steps in order to obtain 4 in the results. The sum of the 
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squared deviations is divided by the number of events (less 1, because devi- 
ations are from the computed mean rather than from the true mean) in order 
to obtain the variance. The square root of the variance is the standard devi- 
ation 

In the fourth step, the numbers of observed values lying within ranges de- 
lineated by various percentages (30, 50, 90, etc) of the standard deviation 
from the mean were counted. The numbers counted were tabulated and com- 
pared with the numbers expected in each range for a normal distribution. 
This comparison is shown graphically in Figure 2, where a normal cumulative 
frequency curve is compared with that determined by the values counted as 
outlined above. The comparison is so close that it can be concluded at this 
point that the aggregate of the values for a great number of stations widely 
distributed geographically follows a normal distribution. The accepted chi- 
square test of goodness of fit cannot be applied strictly in a case such as 
this, inasmuch as not all of the values are entirely independent of each other. 

In the fifth step, it is determined whether the indicated distribution func- 
tions of individual stations vary characteristically or randomly from normal. 
It is entirely possible that some of the regions are characterized by peculi- 
arities of skew or kurtosis,3 even though the aggregate is not. For this pur- 
pose, the 45 stations having 40 or more years of record were selected for 
special study. The latest 20 years of record and the preceding 20 years of 
record for each station were tabulated separately. Deviations from the 
20-year mean were computed for each value (logarithm). Sums of the squares, 
cubes and fourth powers of the deviations were computed for each of the 
ninety 20-year records. Coefficients of dispersion, skew, and kurtosis were 
determined by the following equations: 


(2) 


Sk (3) 


K = (4) 


Computed values of mean, dispersion coefficient, skew coefficient and kurtosis 
coefficient are listed in Table 3. Coefficients of determination between the 
respective pairs of columns are .89 for dispersion, .01 for skew, and .00 for 
kurtosis. Thus, only one percent of the skew evidenced in one period was 
consistently evidenced in the second period. Kurtosis persistence is zero 

(to two places). 

The results obtained in the above analysis are believed to constitute ade- 
quate proof that in the great majority of distributions of annual maximum 
mean-daily flow, the assumption of logarithmic normality is justified for all 
practical uses. 


Graphical Presentation 


For hydrologic design purposes and for simplicity in many applications, 
frequency distributions are most conveniently represented by cumulative 
frequency curves. For this purpose, recorded events are numbered in the 
order of descending magnitude and each magnitude is plotted against the 
relative frequency with which it is exceeded. 

There is a lack of general agreement as to the proper method of plotting 
the points. One line of thought is that the largest event in 20 years has been 
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equalled or exceeded once in 20 years and should therefore be assigned an 
exceedence interval of 20 years. This reasoning is not rigorous, because it 
can similarly be reasoned that all events have been equal to or smaller than 
the largest recorded, giving it an exceedence interval of infinity. 

The second line of reasoning is that the largest event is representative of 
largest events of other similar periods, i.e., that half of the 20-year periods 
at a location will have maximum events smaller than the one in question. 
This appears to be quite reasonable (accepting the inevitable assumption that 
the record is a random selection) and results in the conclusion that the largest 
event of a 20-year record will be exceeded in only half of all 20-year periods 
or somewhat oftener than once in 40 years. (Some 20-year periods may have 
2 or more larger events.) The plotting positions given in Table 4 are based 
on this reasoning. 

An example of graphical construction of a cumulative-frequency curve of 
annual flows is illustrated in Table 2 (b) and Figure 3. The logarithmic- 
probability grid was designed by Hazen such that a cumulative - frequency 
curve of a logarithmic-normal distribution would plot as a straight line. The 
line of best fit determined in Table 2 (a) is obtained simply by plotting the 
antilogarithm of the computed mean at 50 percent and the antilogarithm of 
the sum of the mean and standard deviation of 15.87 percent. The line de- 
termined by these two points is theoretically the best possible estimate of 
the true relation, judging by the plotted data only. 

Of course, it is not necessary to plot such a curve simply to establish a 
frequency value. Values computed analytically are given in Table 2(c) and 
should agree exactly with those read from Figure 3. 


Accuracy of Individual Determinations 


The determination that flood flows can be expressed as a normal distribu- 
tion function makes possible the application to hydrology of reliability criteria 
developed by Student? and others. Student showed that the distribution of 
errors of calculated means from true means follows a particular law repre- 
sented by his “t-distribution”. Similarly, he presented criteria for describ- 
ing the distribution of calculated values of dispersion relative to the true 
value (developed by Pearson and known as the “chi-square distribution”). 

Inasmuch as the true value of a statistic of a continuous variable never 
can be determined exactly, it never will be known how much a calculated 
value differs from the true value. However, limits on either side of the true 
value can be established within which any desired percentage of calculated 
values will lie. This percentage will be determined as a matter of policy in 
any particular application, but for convenience of discussion, a value of 90 
percent will be used in the following illustrations. Five percent of the values 
will lie outside of each of the limits. 

The error of a calculated mean is expressed as the following function of 
the calculated standard deviation: 


(5) 


in which, M = calculated mean 
true mean 
t = function developed by Student 


The value of (M - uw) calculated by use of Equation (5) is exceeded with the 
same frequency that the adopted value of t is shown by an appropriate table? 
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to be exceeded (disregarding sign). Positive values of (M-,) or (u-M) are 
exceeded with half that frequency. 

Similarly, the error of a calculated standard deviation is expressed as the 
following function of the calculated standard deviation: 


(8) 


in which, o = true standard deviation 
X? = function developed by Pearson 


The value of (S - c) calculated by use of Equation (6) is exceeded with the 
same frequency that the value of X? is shown by an appropriate table? not to 
be exceeded. 

Knowing the error of the mean value and the error of the standard devia- 
tion, it is possible to compute the error of a value at any number (x) of 
standard deviations from the mean by the equation: 


6, = (xS - xa)? (7) 


Suppose it is desired to determine the amount to be added to the one-per- 
cent flood calculated from a 20-year record so that the resulting value has 
only one chance in 20 (5 percent) of being smaller than the true one-percent 
value. From published tables, the values of t and xX? correspondiag to .05 
for a series of 20 events are 1.729 and 10.117, respectively. From equation 
(5), (M-u) = -.387 5. From Equation (5), = -.3705. A table? of 
normal-curve values will show that for the one-percent value, x = 2.326. 
From equation (7), 5x = .944 5. A common value of § for logarithms of flood 
flows is .2, in which case 6x would be .189. Addition of this logarithm is 
equivalent to addition of 55 percent to the arithmetic equivalent of the cal- 
culated one-percent value. 

By similar calculations, the errors of a one-percent flood exceeded with 
various likelihoods when using a 5-year record were calculated and plotted 
as Figure 4. To illustrate the applicability of this relationship, records of 
the 10 stations having 50 or more years of record were divided into 5-year 
periods. Inasmuch as the error of the one-percent value obtained using the 
entire record in each case is small compared with the usual errors obtained 
using 5-year records, the use of the long-period determination as the true 
value should not introduce great error for the purpose of this test. Accord- 
ingly, one-percent values were determined using each of the 123 five-year 
periods, and errors measured from the assumed true values were tabulated 
in terms of the standard deviation calculated from the 5-year record. These 
results were arranged in the order of descending magnitude and plotted as a 
cumulative frequency curve in Figure 4. The comparison with the theoretical 
curve appears to be satisfactory. 

A less direct but more rigorous test (inasmuch as it does not involve the 
assumption of a true value) is illustrated in Figure 5. The data of Table 2 
were used to calculate one-percent floods for each 20-year record. The 
difference of each pair of determinations was divided by the average standard 
deviation of the two periods. The resulting values were plotted as a cumula- 
tive frequency curve and compared to the theoretical distribution. The error 
theory appears to apply perfectly. 

From the information presented thus far, it is believed that the theoretical 
dependability criteria which are accepted generally in statistical work are 
fully applicable in the treatment of runoff data. As a practical illustration 
of the general magnitudes of errors involved, Figure 5 shows the 1-in-20 
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limits of error for various record lengths and frequency determinations 
where the calculated standard deviation is .20. The lower curves delineate 
the limits below which the true values will fall with a frequency of once in 20. 
The upper curves delineate the limits above which the true values will fall 
with a frequency of once in 20. Between the corresponding two curves, 9 out 
of 10 of the true values will fall. It can be seen, for example, that when the 
standard deviation calculated from a 50-year record is .20, the true 50- 
percent (mean) value has 19 chances in 20 of being smaller than 1.11 times 
the calculated value, and the true 1-percent value (100 year flood) has 19 
chances in 20 of being smaller than 1.27 (3.72/2.92) times the calculated value. 

Figure 7 is presented for the conveniences of those who are accustomed 
to thinking in terms of probable error. In this case, there is one chance in 
two of the true value lying between corresponding upper and lower curves, 
and one chance in 4 that it will lie above the upper curve. 

Calculation of error-limit curves can be facilitated by the use of Table 5, 
in which values of probable error and 1-in-20 error in terms of the standard 
deviation are given for various exceedence frequencies and lengths of record. 


Need for Region-wide Analyses 


It takes only a glance at Figure § to show that results obtainable from data 
of an individual station are wholly inadequate to the needs of flood-control 
design. This fact has of course been recognized for years. Enveloping curves 
have been developed, intricate hydrologic analyses have been devised, and 
extensive meteorological studies have been carried out for the purpose of 
combining information from many areas and regions for the solution of any 
specific problem. It is necessary to analyze flood-frequency characteristics 
on a region-wide or country-wide basis before adequately precise and reliable 
determinations are possible. 

This problem should be much simpler than the problems that have already 
been attacked by hydrologists, inasmuch as all of the information required 
for the construction of a cumulative frequency curve is summarized in two 
values, i.e. the mean and the standard deviation. These values must be cor- 
related with drainage area characteristics, geographic location, meteoro- 
logical characteristics, etc. The results of such an analysis should make 
possible increasingly accurate determinations for all streams, even for those 
where good records are available. 

Some valuable discussion on the factors to be considered in a regional 
analysis is included in the paper and subsequent discussions on ‘Stream Flow 
Variability” by E. W. Lane and Kai Lei. Error criteria based on Equations 
(5) and (6) are presented graphically in Figure 8 for convenience in making 
regional studies. 


Additional Reliability Considerations 


The analysis developed so far deals with determinations and accuracy 
based on various theories of sampling. Allowances have not been made for 
errors of stream-flow measurements and possible inadequacies of sampling. 

Systematic errors in measurement will generally result in increased 
error of the calculated mean, but will not always affect the accuracy of the 
standard deviation. Their effects should be largely eliminated in a proper 
regionalized study. 

Erratic errors in measurement will affect the accuracy of both of the 
calculated values. However, the effects are much smaller than it might seem. 
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As an example, suppose the largest flood of the 2) shown in Table 4 were 
actually 25 percent larger than was measured. Correcting for this would 
change the mean from 4.944 to 4.949 and the standard deviation from .205 to 
-215, which would increase the calculated one-percent (100-yr) flood by a 
factor of .07. 

The thesis of normality, as applied here, involves the assumption that the 
effects of factors capable of producing a large event are adequately experi- 
enced in the record. In some cases, this may not be true. Dam failures, 
earth slides, and debris jams constitute potential] hazards which are not often 
indicated in the record. As a more common example, Colorado River at 
Grand Canyon experiences an annual flood from snow melt, which registers 
the maximum mean-daily flow almost every year. Yet, summer floods of 
considerable magnitude have been experienced. Their threat is in their rarity, 
and when estimating frequencies on the order of one percent or less, the pos- 
sible effects of such floods should be considered. These conditions should 
result in a positive coefficient of skew in all cases where their effect is ap- 
preciable, which fact, coupled with the observation that there is no predomi- 
nance of positive-skew cases, indicates that the consideration will be re- 
stricted to a relatively few areas. If such cases are found to exist, it will be 
necessary to estimate the effect by region-wide analysis, in which all of the 
extreme floods observed in the region of positive skew are compared with the 
number of extreme floods expected for the case of normality. This type of 
analysis would be advisable in any event as a test of normality. 

Recorded data for most rivers are supplemented by historical data and 
high-water observations for great floods. In the determination of flood- 
frequency relationships for a given location, it is essential that all available 
pertinent information be utilized. This may be done by statistical means or 
as a matter of judgment. In this connection, it is pointed out that statistics 
given in the last two columns of Table 1 are for use in this analysis, but are 
not necessarily the best estimates, inasmuch as historical data were not con- 
sidered. It is doubtful, however, that sketchy historical data would greatly 
modify values based on records of 25 or more years’ length. 


CONC LUSIONS 


1. With rare exceptions, the logarithms of annual maximum mean-daily 
natural flow values are normally distributed. 

2. Effective dependability criteria are readily available for application to 
individual flood-frequency determinations. 

3. At best, flood-frequency determinations based on individual records are 
crude. It will be necessary to make region-wide analyses before precise 
determinations are possible. 
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TABLE I ~~ RIVER STATION DATA 


Computed 
Station Drainage Record Values 


River Location Length 
q-Mi. “Years Tog.j9 Tog-jo 


=) 


Machias R. 
Mattawamkeag R. 
Piscataquis R. 

Little Androscoggin R. 
Pemigewasset R. 


Souhegzan R, 

Br. Westfield 
Housatonic R, 
Schroon R, 

Hoosic R, 


Rahway R. 

Raritan R. 

E. Br. Delaware R, 
Susquehanna 8. 
Chemung R. 


W. Br. Susquehanna R. 
Juniata R. 

Patuxent R, 

Potomac R, 
Rappahannock R, 


Perkiomen Cr. 
James R, 
Appomattox R. 
Roanoke R, 
Dan R, 


Yadkin R, 
Oconee R. 
Chattahoochee R, 
Oostanaula R. 
Tombigbee R. 


Pearl R. 

Allegheny R. 
Blacklick Cr. 
West Fork R. 


Kokosing R. 
Hocking R. 
Reed Cr. 
Greenbrier R, 
Meadow R, 


Whitneyville, Me. 
Mattawamkeag, Me, 
Dover-Foxfort, Me. 
So. Paris, Me. 
Plymouth, N.H. 


Merrimack, N.H. 
Goss Heights, Mass. 
Falls Village, Conn. 
Riverbank, N.Y. 
Eagle Bridge, N.Y. 


Rahway, N.Jd. 
Manville, N.J. 
Fishs Eddy, N.Y. 
Harrisburg, Pa, 
Chemung, N.Y. 


Williamsport, Pa, 
Newport, Pa. 
Burtonsville, Kd. 
Point of Rocks, Md. 
Fredericksburg, Va. 


Frederick, Pa. 
Cartersville, Va. 
Mattoax, Va. 
Roanoke, Va. 

So. Boston, Va. 


Wilkesboro, N.C. 
Dublin, Ga. 
West Point, Ga. 
Resaca, Ga. 
Columbus, Miss. 


Jackson, Miss. 

Red House, N.Y. 
Blacklick, Pa, 
Butcherville, W. Va. 
Rockville, W. Va. 


Millwood, Ohio 
Athens, Ohio 
Grahams Forge, Va. 
Alderson, W. Va. 
Nallen, W. Va. 


3-748 
4.232 
3-871 
3.252 
4.205 


3.492 
3.172 


2117 
2180 
0215 


2196 


457 34 
1,400 
286 A 
76 
622 60 
171 37 | 
3 3% an 
632 34 3.772 
527 39 3.697 162 
510 180 
41 25° 2.881 
R 490 25 4.027 .155 
13 783 34 4.265 
24,100 56 «5.418 147 
15 2,530 43 4.596 .178 
16 5,680 52 4.939 
17 3,350 49 4651 169 
18 127 28 3.164 .219 
19 9,650 50 4.998 
20 1,600 39 4.358 .278 
21 152 29 3.683 .120 
102 6,240 48 4800 .176 
103 729 25 30843 255 
105 2,730 29 
106 493 32 «32958 
108 3,550 50 4611 .210 
109 1,610 50 4.302 .199 
110 4,490 29 4.589 
3,100 27 2552 
201 1,690 43° 4e350 6143 
202 390 42 3.961 .179 
203 181 346674 153 
204 Big Sandy Cr. 200 34 3.813 .182 
| 206 454 25 3.849 231 
207 31 4-091 164 
208 24,7 27 30462 
209 1,360 50 4.415 
210 287 26 32723 
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TABLE I -- RIVER STATION DATA (Conta. ) 


River 


Raccoon Cr. 
Mad R, 

Kentucky R. 
Embarrass R. 
W. Fk. White R. 


Cumberland 


S. Fk. Cumberland R. 


Red R. 
French Broad R, 
Davidson 


S. Fk. Holston R. 
Tuckasegee R. 
Powell R. 
Hiwassee R, 
Sequatchie R. 


Elk R. 
Buffalo R. 
Tennessee R. 
Pike R, 

Fox R. 


Wolf R. 
Milwaukee R. 
Auglaize R. 
Canaseraga Cr. 
St. Regis R. 


Red R. of North 
Red Lake R. 
Pembina R. 
Mississippi 
Crow R. 


Minnesota R, 
Black R, 
LaCrosse R, 
Root R. 
Turkey R. 


Maquoketa R. 
Rock R, 
Sugar R. 
Iowa R, 
Cedar R, 


Station 
Location 


Adamsville, Ohio 
Springfield, Ohio 
Winchester, Ky. 

Ste. Marie, Ill. 
Noblesville, Ind. 


Nashville, Tenn. 
Nevellsville, Ky. 
Adams, Tenn. 
Asheville, N.C. 
Brevard, N.C. 


Bluff City, Tenn. 
Bryson City, N.C. 
Arthur, Tenn. 


Murphy, N.C. 
Whitwell, Tenn. 


Prospect, Tenn. 
Flat Woods, Tenn. 
Chattanooga, Tenn. 
Amberg, Wis. 
Berlin, Wis. 


New London, Wis. 
Milwaukee, iis. 
Defiance, Ohio 

Dansville, N.Y. 


Brasher Center, N.Y. 


Grand Forks, N. Dak. 


Crookston, Minn. 
Neche, N. Dak. 
Clinton, Iowa 
Rockford, Minn. 


Mankato, Minn. 
Neillsville, Wis. 
W. Salem, Wis. 
Houston, Minn. 
Garber, Iowa 


Maquoketa, Iowa 
Como, Ill. 
Brodhead, Wis. 
Marshalltown, Iowa 
Cedar Rapids, Iowa 


Computed 


Drainage Record 


Values 


Years log.jo9 10g-10 


3-806 
3-686 
4.701 
4.024, 
3.979 


52076 
4.626 
4.052 
4.140 
30043 


4.029 
4.034 
4.107 
30933 
3.897 


30985 
50298 
2.998 
32563 


30787 


2157 


| 
Sq... 
2u1 587 28 
212 8532 
213 3,960 37 
214 1,540 33 1263 
215 845 31 2229 
q 
217 12,900 109 
218 1,260 30 
219 678 (26 
221 40 25 e211 
222 816 46 
2h 685 27 2165 
225 406 2219 
227 1,780 28 .170 
228 LA? 26 
229 21,400 61 
301 250 3 
302 1,430 1g 
303 2,240 50 .194 ' 
304, 661 32 34634 
306 2,330 26 4e305 
307 153 28 3.376 .186 
308 616 33 30840 .126 
401 30,100 49 4.040.328 
4,02 5,280 45 3.689 .310 
4,03 3,080 39 3.029 .413 
88,600 73 5.113 .153 
1,05 23520 25 1,328 1463 
1,06 14,900 
407 756 37 4.030 
1,08 206 
4,09 1,270 25 32908 .260 
410 1,530 25 4.082 .180 
411 1,550 33 4.058 222 | 
412 8,700 32 4.366 0232 
413 529.323.520.305 
1,530 25 36775 
415 6,640 43 
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River 


Skunk R, 

Des Moines R, 
Kankakee R. 
Fox R. 
Mackinaw R. 


Sangamon R. 
Big Muddy R. 
St. Croix R. 
Madison 
Lodge Cre 


Republican R. 
Kansas R, 
Solomon R, 
Big Blue R. 
Grand R. 


Medicine Cr. 
South Grand R. 
Yellowstone 
Meramec R. 

St. Francis R. 


James R, 
Eleven Point R. 
Walnut R, 
Neosho R, 

Petit Jean Cr. 


Sulphur R. 
Sabine R. 
Neches R, 
Trinity R. 
Brazos R. 


Little R. 
Colorado R. 
N. Llano R, 
Guadalupe R. 
Conejos R. 


Blue R, 

Elk R. 

Ashley Cr. 
Huntington Cr, 
Tonto Cr, 


Station 
Location 


Augusta, Iowa 
Keosauqua, Iowa 
Momence, Ill. 
Dayton, Til. 
Green Valley, Il. 


Monticello, I1l. 
Plumfield, Ill. 
Danoury, Wis. 

W. Yellowstone, Lont. 
International Bdy. 


Clay Center, Kan. 
Topeka, Kan. 
Niles, Kan. 
Randolph, Kan. 
Gallatin, Mo. 


Galt, Mo. 
Brownington, lo. 
Glendive, Mont. 
Eureka, Mo. 
Patterson, Mo, 


Galena, Mo. 
Bardley, Mo. 
Winfield, Kan. 
Iola, Kan. 
Danville, Ark. 


Darden, Tex. 
Logansport, La. 
Rockland, Tex. 
Oakwood, Tex. 
Waco, Tex. 


Cameron, Tex, 
Ballinger, Tex, 
Junction, Tex. 

N. Braunfels, Tex. 
Mogate, Colo, 


Dillon, Colo. 
Clark, Colo. 
Vernal, Utah 
Huntington, Utah 
Roosevelt, Ariz. 


TABLE I -- RIVER STATION DATA (Contd,) 


Drainaze Record 


Area_ Length 


Sq.M. 


13,900 39 
2,340 32 
2,570 32 
1,100 25 


550 

753-36 
1,550 

419 

797 


24.4600 
56,700 
6,770 
9,100 
2,250 


Years 


225 
1,660 
66,100 
3,788 
956 


987 


Computed 
Values 


log.19 10g-j0 


No. |_| | | 
416 
417 210 
418 30739 163 
419 30940 2252 
420 30833 32K 
421 30639 304 
422 30837 
42h 
502 3.054 123 
508 320050 348 
519 40 4.078 
520 LO &.5€7 313 
522 35 3.776 .346 
523 28 4.262 «310 
524 25 .303 
525 25° 3.518 3 
526 25 4.079 .358 
527 4h 
601 25 4505 284 
602 25 2326 
603 |_| 25 320 
604 793 25 32756 467 
606 1,890 25 e218 .341 
607 3,800 48 
609 760 29 
610 2,750 23 313 
651 4,860 39 319 
652 3,540 4O 4.146 
653 12,800 34h 4.509 371 
654 19,300 42 e666) 6286 
655 7,030 30 
565 5,340 38 
657 914 31 «630401 «w. 7555 
658 1,666 31 3.810 .532 
663 282 4h 30385) 
704, 129 7 28s 
709 206 27 119 
710 29 2.879 4.196 
711 188 22, 
716 678 34 3547 499 
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TABLE I — RIVER STATION DATA (Contd 


River 


Blacksmith Fk. 
Jeber R. 
Spanish Fk. 

W. ‘ialker R. 
Chewaucan I, 


Silvies 
Santa Ysabel Cr. 
Santa Ana R, 
Waterman Cn. Cr. 
San Gabriel R. 


L. Santa Anita Cr. 
Arroyo Seco 

Kern R. 

Tule R. 

Kings R. 


Ne Fx. Kings.R. 
Tenaya Cr. 

N. Fk. Feather 
Oregon Cr. 
Quinault Re 


S. Fk. Skykomish R, 
Colwabia R, 
Kootenai 

St. Joe a. 

Boise R. 


Salmon 2. 
Grande Ronda 
Clearwater i. 
Columbia 8, 
Crooked R, 


little Sandy R. 
Willamette R. 
McKenzie R, 
Clackamas R,. 
Cowlitz R. 


Big Butte Cr. 


Station Drainage Record 
Location Area_ Lengi:h 
Sq-Mi. Yeai's 


Hyrum, Utah 260 
Oakley, Utah 

Thistle, Utsh 

Coleville, Calif. 

Paisley, Ore. 


Burns, Ore. 

Mesa Grande, Calif. 
Mentone, Calif. 
Arrowhead Springs, Calif. 
Azusa, Calif, 


Sierra Madre, Calif. 
Soledad, Calif. 
Bakersfield, Calif. 
Porterville, Calif. 
Piedra, Calif. 


Cliff Camp, Calif. 
Yosemite, Calif. 
Bidwell Bar, Calif. 
San Juan, Calif. 
Quinault Lake, Wash. 


Index, iashe 
Birchbank, B.C, 

lil DY» Mont. 

Calder, Idaho 

Twin Springs, Idaho 


Whitebird, Idaho 
La Grande, Ore, 
Kamaiah, Idaho 
The Dalles, Ore. 
Culver, Ore. 


Bull Run, Ore. 
Albany, Ore. 
McKenzie Vridce, Ore. 
Cazadero, Ore. 
Packwood, 


Siletz, Ore. 
Elkton, Ore, 
Powers, Ore. 
Butte Falls, Ore. 


Computed 
Values 


10g +16 


: 
752 33 2.570 
30 2.685 
757 3.197 187 
760 3h 2.813 
761 93h 31 30954 0349 
801 58 22 2.613 2585 
803 195 50 2.780 .642 
804, 5 30 1.581 .426 
£06 211 38 30296 2668 
607 2 30 2601 
808 241 45 30580 390 
809 2,420 49 36571 0290 
210 266 30156 
1,690 51 217 
812 17h 25 3ehl5 
816 1,350 35 
817 35 35 290 
901 264 31 0173 
902 355 37 218 
904, 10,200 36 4he751 157 
906 1,080 25 4ell6 
934 830 36 30759 2160 
935 12,600 3h e733 e148 
936 678 36 3431 
937 4,850 36 0136 
961 237,000 77 56782 .127 
963 39830 29 3.619 172 
964 22 30129 185 
965 4y8hO 999 206 
966 365 35 30718 .198 
967 665 37 «e252 227 
969 Siletz R. 202 27 =e 250 2150 | 
970 Umpqua R. 3,680 41 64846 
971 S. Fk. Coquille R. 169 36990 «6.155 
973 Fk. 131 29 22807 
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TABLE 2 — VALUES FROM TYPICAL 20-YEAR RECORD 


(Willamette River at Albany, Oregon) 


Table 2(a) 


Water Maximum M.O.F. 


Year log. 
10 
c.f,s 


Exceedence interval 
Exceedence frequency 
3S mltiples above mean 
Diff. from mean (log.) 
Mean (log.) 

Computed flow (log.) 


Computed flow (cfs) 


489 «7495 


488.9517 
0.797E 
0.04199 
0.205 


Taole 2(b) 


Antilog. = 141,000 
Plot at 15.9%. 


M -S = 4.739. 
Antilog. = 54,800. 
Plot at 84.1%. 


Flows Position 
efi. 
1927 191,000 5.281 210, 000 
1928 7,000 191,000 8.3 
1929 50,800 he 706 1€5,000 1342 
1930 84,100 e925 134,000 18.1 
1931 108,000 5.033 124,000 23.0 
1932 134,000 50127 120,000 2709 
1933 90,000 e954 108, 000 32.8 
1934 70,200 42846 92,600 3707 
1935 72,600 4.6861 90,000 
1936 120,000 52079 8,100 4705 
1937 124, 000 52093 83,400 5205 
1938 92,600 4.967 744,000 570k 
1939 58,900 40770 724600 62.3 
1940 58, 300 42766 71,200 67.2 
1941 40,100 1.603 70,200 72.1 
1942 £3,400 4e921 58,900 770 
1943 210,000 50322 58, 300 81.9 
19h y4,00 50,800 $6.8 
1945 71,200 Ady 400 9167 
1946 185,000 50267 40,100 96.6 
=X 98.889 
M 
2-yr | S-yr | lyr | 25yr | 50-yr | 100-yr 
0,50 0.20 0,10 0.04 0.02 0,01 
0,842 1,282 1.751 2.054 2.326 
96173 | 0.263] 0.359] 0.421] 0.477 
| 50117} 50207] 56303 | 54365] 5.421 
37,900 | 133,000 | 161,000 201,000 232,000 | 264,000 
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TABLE 3 -- COMPARISONS OF STATISTICS FOR SUCCESSIVE PERIODS 
Latest 20 Years Preceding 20 Year's 
SK 


s SK M 3 


0125 258k 0463 
0789 0485 
0222 1,007 018 
0196 -.960 


0195 1.549 5 309 
2289 2156 -.109 
0225 -.194 -.100 
0259 024 -2950 


-.436 
2652 
2912 

-.051 


-.173 
2322 
-2705 
2185 


-.416 


NENNN 
w Ww Oo 


ee 


STA ee K ‘ 
2 4.219 3.02 
3 3.216 2.47 
5 4.200 2.25 
5.397 2.17 
15 42620 2.87 
16 40909 3054 
17 42626 265k 
19 5.030 3.16 
102 1.89 
104 30734 2.78 | 
108 588 2.34 228 01300 291 
109 4.367.212 3.45 e332 «192 
201 4.350 131 2.40 4e350 151 1.64 
202 3.941 193 3051 30975 «162 20502 2.67 
209 4.421 197 2686 196 -.129 2.02 
217 5.072 2.10 5.058 .090 2&1 
222 3.975 17h 3-35 4.051 -.105 1.93 
223 3.971 2.60 4.036 e210  2.8h 
225 3.861 1.84 3.960 2th 
29 5.208 .127 3.08 123 22s 
302 36549 «200 -.552 2.81 3.600 0155 2091 1.62 
303 3-831 192 -.933 350 3.803 .186 0037s 2 11 
401 309k 14006 -.630 2.70 4.069 969 
4,02 36605 2305 9h 3.722 29K -.615 
4,04, 5.067 0153-16220 he 5.086 .137 
1,06 30954 339 -.967 3.61 4.067 390 -.401 
415 333 -.484 3.05 188 
519 4.092 3,56 4.063.346 
520 4.586 358 -.248 1.87 252 412 
527 4.777 -.681 3054 4.908 .168 -1.220 
607 he 32h -.687 2,65 4.280 
652 4.186 323 312 1.84 4.106 2389 2593 
753 3166 253 3.372 «139 .019 
757 36143 185 e727 3.250 -.155 
803 2.751 0313 2643 2.832 .658 -.122 
809 3.562 2319 3.627 0293 
810 32125 -.056 1.98 2.330 
811 4.079 2229 362 2.96 4.103 2192 -.693 
961 5.720 e141 -.602 2-11 5.783 2090 -.638 
965 he 0287 2.29 5.016 .198 -.193 
970 4.869 176 e479 4.846 -.167 
808 3. 576 2375 -1. 34,0 4.60 3.620 -.204 
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percent) = 


: 


TABLE 


Number of Events 


50.0 46.5 43.5 


9.0 
7.0 
0) 
1 
el 
ol 
1 
1 
2 
2 
2 
2 
3 
3 
3 
3 
3 
3 
4 


bel 506 502 hed 
14.9 13.7 12.7 11.8 11.0 
25.9 23.7 21.8 20.1 18.7 17.5 
5.6 29-8 27.6 25-7 2h.0 
5.2 hle2 37-9 3501 32.6 30.5 
50.0 46.0 42-5 39.6 37.0 

13. 

i. 


0.0 
o2 
0.3 
8.3 
o4, 
5 
5 
6 
6 
7 


6 
16 


30. 
28. 
26. 


4 
20 
1 
05 
26 3 
7 
8 3 
9 
0 
0 
1 


3 
3 
356 
3 
31. 


9 
7 
5 
30 
1 
0 
8. 


AN AONTH 


SSESSIRRRAAS 


2869 294 


2502 2507 26. 
ol 21.5 21.9 22.3 22.8 23.2 


19.3 1967 20,0 20.4 20.8 21, 
o2 17.5 17.8 18.2 18.5 18.9 19. 
e4 1567 16.0 16.3 16.6 16.9 17. 


2303 2308 2he2 Zhe 
06 1369 Woh 14.7 15.0 15. 


24 
22 
21 


e9 12.1 12.3 12.5 12.8 13.1 13. 

10,1 10.3 10.5 10.7 10.9 11.1 11. 

8.3 8.5 8.6 88 9.0 91 

604 6.6 667 668 6.9 7.0 702 7. 

500 500 Sel 502 5e 

209 300 3.0 Fel Bel 302 3.2 3.3 36 
162 1.2 1.2 163 1,3 2.3 14 


CN ON ON OY R 


49-1 
4 
3 
0.6 
9 
ok 
20 
6 
9 
02 
20 
6 
302 
08 
20 


48 47 


Number of Events 


57 


58 


/n) 


100 (1- .5 2 


Other plotting positions interpolated linearly to 50 percent 
Plotting positions are symmetrical about 50 percent 


NOTE: Plotting position for m= 1 obtained from formla, 


i 
TONS 
1 2 3 4 5 6 7 8 i 22: m 
50.0 29.3 20.7 15.9 13.0 10.9 94 8.3 7d 1 
50.0 38.6 31.5 26.5 22.9 20.2 18.1 2 
30 = (49. 50.0 42.2 36.5 32.1 28.7 3 
29: &7- 3 50,0 50.0 44.0 39 4 
28 45e 6 48.3 1 50.0 50 5 
27s ithe 9 46.5 3 48.2 6 
6 3 44.8 6 46.4 7 
25 40, 6 43,0 8 44.6 8 
2, 39. 9 41.3 O 42.8 9 
23 376 2 5 2 41.0 10 
5 38.5 39.2 41 nu 
8 ol 36.7 
20 32. 2 03 3409 3506 13 
19 le 5 «6 33.2 33.8 
18 29. 8 oS 31-4 31.9 15 
a7 1 el 29.6 30.1 
4 16 26. 4 3 27.8 ie 
15 06 26.1 28.5 27.0 27.5 2 
wy 22. 1 09 24.3 
13 4 el 22 
hk 19 7 e4 20 
0 18.6 19 
» 3 16.9 17 
9 14. 6 15.2 15 
12. 13.4 13 
7 3 
6 
5 Te 9 
4 6. 2 
3 4. 5 
2 2e 8 
1 1. 2 
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F 


(percent ) 


Number of Events 


OS 


e2 28.1 27.1 26.2 25.3 
09 31.7 30.6 2926 28.6 
07 3504 Bhel 3320 3129 
05 3900 3707 3604 352 


an< 


an 35 


$2999 424978 
AANRARYS 


ASI 


OO 


dds 


HOH 


ss 


339 3999299799 


$9 


SOAS id 


so a 


ONO Ad 
NAANNN 


33.32 


45 


Number of Events 


Other plotting positions interpolated linearly to 50 percent 
Plotting positions are symmetrical about 50 percent 


NOTS: Plotting position for m* 1 obtained from formla, 100 (1 - .5 1/n) 


(Contd. ) 
4 3e 2.6 265 204 2-3 
10.3 9 3. 604 6.2 6.0 568 526 
30 16.4 15 13. 10.2 968 905 869 
29 22.5 21 19. 14.0 13.5 13.0 12,6 12,2 
27 = 32 296 21.6 20.8 20,0 19.4 18.7 
26 40.8 She 25eh Zhe 2306 22.8 22,0 
25 46.9 396 29 
24 Ae 32 
23 50,0 50. 36 
22 «47.8 48.9 50.0 
21 45.6 46.6 47.7 3 42.7 41.2 39.8 38.5 
20 45-4 8 1 46.3 43.2 41.8 
19° 42.1 43.1 6.3 50.0 48.2 46.6 45.1 
18 39.0 39.9 40.8 308 50.0 48.4 
17 6 38.5 3 
16 34.6 2 
15 32h 9 4, 
30.2 8 6 9 
13 6 3 
12 25.8 3 0 0 
11 23.6 1 6 5 
10 21.4 8 3 0 
9 Bee 6 0 5 
8 17.0 3 7 0 
7 
6 222.5 
5 10.3 6 8 6 
4 8.1 3 5 
3 509 2 
1 1.5 6 6 
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ST 


Averoge volves of 
159 stations (5726 FIGURE 2 


Station-Yeors) GOODNESS OF FIT 


(STANDARD DEVIATIONS FROM MEAN) 


LOGARITHM OF ANNUAL MAXIMUM MEAN-DAILY FLOW 


20 30 4 50 60 70 8 
EXCEEDENCE FREQUENCY -% 


FIGURE 3 


EXAMPLE 
CUMULATIVE FREQUENCY 


ANNUAL MAXIMUM MEAN-DAILY FLOW IN CFS 


20 30 40 SO 60 70 80 
EXCEEDENCE FREQUENCY -% 


| ! 
3 
| 
| | 
aN 
| 
| | | | | |_| | 
oo os 2 5 10 | 95 995 
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COMPARISON OF OBSERVED 
AND THEORETICAL ERRORS 


S~YEAR RECORDS 
ONE-PERCENT FLOODS 


Errors of colculated vawes 
from assumed true volves 


ERROR (STANDARD DEVIATION MULTIPLES) 


20 4 6 70 80 99 6 98 9 995 999 


EXCEEDENCE FREQUENCY—% 


FIGURE 5 


COMPARISON OF OBSERVED 
AND THEORETICAL ERRORS 


20-YEAR RECORDS 
ONE-PERCENT FLOODS 


ERROR (STANDARD DEVIATION MULTIPLES) 


EXCEEDENCE FREQUENCY —% 


Theoretical! distribution 
ary 
| | } 
| | | 
| 
|| | | 
| 
F 
| | 
6 Theoretical distribution 
of differences 
12 T 
Differences between indepen- 
dently colculoted values 
= | 
{ 2 | | | | | 
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FIGURE 6 


THEORETICAL 


value | ) RROR LIMITS 
| | ¥ 
! 


Lumits exceeded by true 
__ value in 95 percent of 
cases 


FLOW (RATIO TO 5O-PERCENT VALUE) 


60 
EXCEEDENCE FREQUENCY -% 


Limits exceeded by true valve in 
25 percent of oll case: 


ccleuloted volve_| FIGURE 7 


($= 0.2) THEORETICAL 
ERROR LIMITS 


Limits exceeded by 
true in 75 per- 


—— cent of all coses - 


FLOW (RATIO TO 5O-PERCENT VALUE) 


| } 


| 
20 30 40 50 60 70 60 ” 
EXCEEDENCE FREQUENCY—% 


NET 5 percent of aii cases 
_ 
| | | 
| | | 
oO. os 1 2 5 90 95 98 995 999 9999 
| 
| 
f 
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